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Abstract 

The Kardar-Parisi-Zhang (KPZ) equation in 1+1 dimension dynamically de- 
velops sharply connected valley structures within which the height derivative 
is not continuous. There are two different regimes before and after creation 
of the sharp valleys. We develop a statistical theory for the KPZ equation in 
1+1 dimension driven with a random forcing which is white in time and Gaus- 
sian correlated in space. A master equation is derived for the joint probability 
density function of height difference and height gradient P{h — h, dxh, t) when 
the forcing correlation length is much smaller than the system size and much 
bigger than the typical sharp valley width. In the time scales before the cre- 
ation of the sharp valleys we find the exact generating function of h — h and 
d^h. Then we express the time scale when the sharp valleys develop, in terms 
of the forcing characteristics. In the stationary state, when the sharp valleys 
are fully developed, finite size corrections to the scaling laws of the structure 
functions {{h - h^ [d^h)'^) are also obtained. PACS: 05.70.Ln,68.35.Fx. 
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I. INTRODUCTION 



There has been a great deal of recent work on the formation, growth and geometry of 
interfaces [1-5] . The dynamics of interfaces has turned out to be one of the most fascinating 
and at the same time challenging topics in theoretical non-equilibrium physics. There has 
been two principal approaches for theoretical analysis of such problems. The first is based 
on computer simulations of discrete models and often provides useful links between analytic 
theory and experiments. The second approach aims to describe the dynamical process by a 
stochastic differential equation. This procedure neglects the short length-scale details but 
provides a coarse-grained description of the interface (that is suitable for characterising the 
asymptotic scaling behaviour). Theoretical modelling of growth processes started with the 
work by Edwards and Wilkinson (EW) [6]. They suggested that one might describe the 
dynamics of the height fluctuations by a simple linear stochastic equation. Kardar-Parisi- 
Zhang (KPZ) [7], realized that there is a relevant term proportional to the square of the 
height gradient which represents a correction for lateral growth. Indeed the KPZ equation 
is a prototype model for a system in which the interface growth is subjected to a random 
external flux of particles. The randomness is described by an annealed random noise which 
mimics the random adsorption of molecules onto a surface. In the KPZ model (e.g in the 
1-1-1 dimension), the surface height h{x, t) on the top of location x of 1-dimensional substrate 
satisfles a stochastic random equation, 

§-^id^hr = udih+fix,t), (1) 

where a > and / is a zero-mean, statistically homogeneous, white in time and Gaussian 
process with covariance 

{fix, t)f{x\ t')) = 2Do6{t - t')D{x - x'), (2) 

where, 

D{x-x') = ^exp(-^^^), (3) 
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and a is the variance of D{x — x'). Typically the correlation of forcing is considered as 
delta function for mimicking the short range correlation. We regularise the delta function 
correlation by a Gaussian function. When the variance a is much less than the system size 
we would expect that the model would represent a short-range correlated forcing. So we 
would stress that our calculations are done for finite a ^ L, where L is the system size. 
The average force on the interface is unimportant and may be removed from the equation 
of motion. Every term in the eq.(|I|) involves a specific physical phenomenon contributing 
to the surface evolution. The parameters z/, a and Dq ( and a) are describing the surface 
diffusive relaxation, non-linear lateral growth and the effective noise strength, respectively. 

We consider a substrate of size L and define the mean height of growing film and its 
roughness w by, 

1 /•^/2 

h{L,t) = — dxh{x,t), (4) 

L J-L/2 

^L,t) = (((/^-/^)^))V^ (5) 

where (■ ■ ■) denotes an averaging over different realizations of the noise (samples). Start- 
ing from a flat interface (one of the possible initial conditions), it was conjectured by Family 
and Vicsek [8] that a scaling of space by factor b and of time by a factor {z is the dynamical 
scaling exponent), re-scales the roughness w by factor b^ as follows, 

w{bL,bH) =b''w{L,t), (6) 

which implies that 

wmt) = L^f{^). (7) 

If for large t and fixed L {jj —>■ oo), w saturates then f{x) const, as x ^ oo. However, 
for fixed large L and 1 << t << , one expects that correlations of the height fluctuation 
are set up only within a distance t^^^ and thus must be independent of L. This implies that 
for X << 1, /(x) ~ with = f • Thus dynamic scaling postulates that [9], 
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The roughness exponent x ^'^'^ the dynamic exponent z characterise the self-affine geometry 
of the surface and its dynamics, respectively. Several time-regimes can be distinguished in 
the time evolution of the surface roughness. They can be summarised as follows: for very 
early times , the noise term dominates since its contribution to the equation grows as the 
square root of time. In this time- regime, the surface roughness grows as w{i) ~ i^/^. For 
intermediate times, the linear term has the main contribution. The linear case {a = 0) 
is the Edwards- Wilkinson model for which one can easily find that the surface roughness 
behaves as w{t) ~ t^°, where the value of Po depends on dimension of substrate(/3o = 
^ for d-dimensional surface). For later times, the contribution of the relevant non-linear 
term becomes a dominant one and the surface roughness growth is characterised by the 
behaviour w{t) ~ . For very late times and finite substrate- length L, the roughness 
saturates to the value ^(i — > oo, L) ~ L^. Of course, in an experiment or in a numerical 
simulation the transition between the different regimes is not sharp and different crossover 
behaviours can be observed. Galilean invariance implies the relation % -|- 2; = 2 independent 
of dimension [10-11]. It means that there is only one independent exponent in the KPZ 
dynamics. In the one-dimensional substrates a fluctuation-dissipation theorem yields exactly 
z — ^, X — ^:P— ^ [12]. In contrast to 1-dimension, the case d>2 can be only attacked by 
approximative field-theoretic perturbative expansions [13-16]. It is well-known that effective 
coupling constant for the KPZ equation is g = ^"^3^" . Phase diagram information extracted 
from the renormalisation group flow indicates that d — 2 plays the role of a lower critical 
dimension. For d < 2, the Gaussian fixed point {a — 0) is infrared-unstable, and there is 
a crossover to the stable strong coupling fixed point. For d > 2, a. third fixed point exists, 
which represents the roughening transition. It is unstable and appears between the Gaussian 
and strong coupling fixed points which are now both stable. Only the critical indices of the 
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strong-coupling regime (g^ — > oo or i/ — > 0) are known in 1+1 dimensions and their values in 
higher dimensions as well as properties of the roughening transition have been known only 
numerically [17-23], and the various approximation schemes [24-32]. 

The theoretical richness of the KPZ model is partly due to close relationships with other 
areas of statistical physics. It is shown that there is a mapping between the equilibrium 
statistical mechanics of a two dimensional smectic-A liquid crystal onto the non-equilibrium 
dynamics of the (1+1)- dimensional stochastic KPZ equation [33]. It has been shown in [34] 
that, one can map the kinetics of the annihilation process A + S — > with driven diffusion 
onto the (l+l)-dimensional KPZ equation. Also the KPZ equation is closely related to 
the dynamics of a sine-Gordon chain [35], the driven-diffusion equation [36-37], high Tc- 
superconductor [38] and directed paths in the random media [39-52] and charge density 
waves [53], dislocations in disordered sohds [3], formation of large-scale structure in the 
universe [54-57] , Burgers turbulence [58-85,90] and etc. 

As mentioned the main difficulty with the KPZ equation is that it is controlled, in all 
dimensions, by a strong disorder ( or strong coupling) fixed point and efficient tools are 
missing to calculate the exponents and other universal properties e.g. scaling functions, 
amplitudes, etc. Despite the fact that in one dimension, the exponents are known, but 
many properties, including the probability density function (PDF) of the height of a growing 
interface have been so far measured only in numerical simulations. Recently, it is shown that 
for one particular model of the KPZ class, the asymmetric exclusion process (ASEP), the 
whole distribution of displacement of particles could be calculated for a finite geometry by 
Bethe ansatz [86-87]. As discussed in [86-87], the PDF of displacement of particles is related 
to the PDF of height differences in the growth process. It is proved in [86-87] that the PDF 
of height differences has the following asymptotic: 
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where y — and w is the variance of the height fluctuation about its average. This result 
is in agreement with the numerical experiments [88]. 

In this paper we are interested in the statistical properties of the KPZ equation in 
the strong coupling limit {v 0). The limit is singular, i.e. through which the surface 
develops sharp valleys. So starting with a flat surface after a flnite time scale, tc, the sharp 
valley singularities are dynamically developed. In the singular points (sharp valleys) spatial 
derivative of the h[x, t) is not continuous. Hence the limit of i/ ^ is not singular for t < tc, 
and we can ignore the diffusion term while after developing the singularities the diffusion 
term has finite contribution in the PDF of height fluctuations. Inspired by the methods 
proposed, recently in the works done by Weinan E and Vanden Eijnden [73], we develop a 
statistical method to describe the moments of height difference and height gradient of height 
field h{x,t). We derive a master equation for joint probability density function (PDF) of 
the height difference {h — h) and height gradient dxh, P{{h — h), dxh) for given g or diffusion 
constant v. We will consider two different time scales in the limit of — > 0. (i) Early 
stages before developing the sharp valley singularities and {ii) established stationary state 
comprising with fully developed sharp valley singularities. In the regime (i) ignoring the 
relaxation term in the equation of the joint PDF when i/ ^ 0, we determine the exact 
generating function of joint moments of height and height gradient fields. Realisability 
condition for the resulting joint PDF sheds light on the time scale of the sharp valley 
formation. In contrary the limit — > is singular in the regime {ii) leading to an unclosed 
term (relaxation term) in the PDF equation. However we show that the unclosed term can be 
expressed in terms of statistics of some quantities defined on the singularities (sharp valleys). 
Identifying each sharp valley in position yo with three quantities, namely the gradient of h 
in the position |/o+, Uq- and its height from the we determine the dynamics of these 
quantities. In the both regimes all the moments, {{h — h)^{dxh)"^) , for given n and m are 
found. In the regime (ii) we will prove that in leading order, when L — > oo, fluctuation 
of the height fleld is not intermittent and succeed to give the analytic form of amplitudes 
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of the all structure functions. Besides, the scaling behaviour and the amplitudes of all the 
correction terms due to the finite size effect are calculated. 

The paper is organised as follows; in section two, we derive the master equation for the 
joint PDF of height differences and height gradients for given diffusion v. We convert the 
height PDF, i.e. P{h — h,t), evolution equation consequently to a Fokker-Planck equation 
for an arbitrary given diffusion constant. In section three we will consider the limit of ^ 
of the master equation in the time scales that there are no any singularities in the surface 
(before developing the sharp valleys). We determine the exact and explicit expression of 
the generating function for the moments {{h — hy'' (d^h)"'') for given n and m. In section 4 
we consider the master equation in the limit i/ — > and consequently when the singularity 
are fully developed. In this regime the relaxation term has finite contribution in the master 
equation. Using the methods introduced in [73], we prove that the unclosed term can be 
written in terms of quantities which are defined on sharp valleys where dxh is discontinuous. 
Also in this section we determine the relation between the density of sharp valleys and the 
forcing variance kxx{0) in detail. In section 5, we derive the moments of height fiuctuation 
in the stationary state and show that the PDF of {h — h) is strongly asymmetric and prove 
that in to leading order the n-th moments of {h — h), i.e. {{h — h)"'), can be written in 
terms of the second order moment of height fiuctuation in a non-intermittent way. We 
determine the amplitudes of the all of moments and show that the amplitudes of moments 
{{h — hY{dxh)'^) can be written in terms of characteristics of singularities. We also derive 
the finite size effect on the moments of height differences and determine the amplitudes of 
all correction terms. In section 6, we derive the PDF of quantities which characterise the 
singularity, and therefore find the equation for their evolution. 
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II. THE MASTER EQUATION FOR HEIGHT DIFFERENCE AND HEIGHT 

GRADIENT 



In this section we consider the 1 + 1— dimensional KPZ equation and derive the master 
equation to describe the joint-PDF of height-difference and height-gradient i.e. P{h~h, d^h), 
for given u and a. It is shown that the equation for the joint PDF is not closed due to 
the linear term ud'^h. The PDF of height difference is related to the joint PDF, P{h — 
h,d^h,t) by the relation P{h -h,t) = - h,d^h,t)d{d^h). We show that P{h - 

h, t) satisfies a Fokker-Planck equation and write down the explicit expression of drift and 
diffusion coefficient D^^^ and D^'^\ It is shown that the drift and the diffusion coefficients 
can be written in terms of the conditional average of {{dxh)'^\h — h). 

We consider a one-dimensional line of length L and a surface of height h{x,t), and 
gradient dxh{x,t) at time t. The 1+1 dimensional KPZ equation governed over height field 
h{x, t) is defined in eq.(|l]) while u{x, t) = —d^h^x, t) is a solution of the the so called Burgers 
equation as, 

ut + auu^ = vu^^ - f^{x, t), (10) 

where the covariance of / is given by eqs.(0) and (|]). To investigate the statistical properties 
of eqs. (p!OD and (0), let us define the generating function Z{X, fi,x,t) as: 

Z{X, fi, x,t) = {exp{—iX{h{x,t) — h) — ifiu{x,t))). (11) 

It follows from eqs.(^) and ([1|) that generating function Z is a solution of the following 
equation, 

z,= i^{t)\z - z,, - x^mz + ^^^±^z, - 
z 

-iafx{ — )^ + fi'^k^^{0)Z - ifiu{u^x exp{-iX{h - h) - ifiu)), (12) 
/i 

where k{x - x') = 2DoD{x - x'),-f{t) = ht, k{0) = ^ and k^^{0) = To derive 

eq.(|12|) we have used the following identities: 



— — % 

{hxx exp {-i\{h{x, t) - h) - i^u{x, t))) = — {Z^ + Xd^Z}, (13) 
(/(x, t) exp {-iX{h{x, t) -Ji) - ifiu{x, t))) = -iXk{0)Z, (14) 

and 

{fx{x, t) exp {-iX{h{x, t) -h) - ifiu{x, t))) = -ifikxx{0)Z, (15) 
where we have used the fact that -0^(0) = 0. This is evident that due to the last term in 



the eq.([l^, it is not closed. Assuming statistical homogeneity {Z^ = 0) we have, 



- ifiZt= jXfiZ - -XfiZ^^ + iX'^fik{0)Z - i^^kxx{0)Z 



—i{vX? + iaX)Z^ — fi^u{uxx exp(—iXh{x, t) — ifiu{x, t))), (16) 

where h{x, t) = h{x, t) — h. Defining P{h, u, t) as the joint probability density function of h 
and u, one can construct P{h, u, t) in terms of generating function Z as, 

P{h,u,t) = [ [ ^^exp{iXh + ifiu)Z{X, fi,t). (17) 
J J 2n Ztt 

It follows from eqs. ([T7|) and (|16D that P{h,u,t) satisfies the following equation, 

- Put= -iP-,^ - ^{u'P)k,u - (^{UP\ - MO)P^^„ + kxxmuuu - HnPh-H 



— exp{iXh + ifiu)^'^{uxx ex.p{—iXh{x, t) — ifiu{x, t))). (18) 

ZTT ZTT 



Now we can rewrite the last term in eq.(|TH[) as following, 

exTp{iXh + ijj,u)iJ^ {uxx{x, t) exp(—iXh{x, t) — ijj,u{x, t))) 



271 2n 



—^{uxx{x, t)6{h{x, t) — h)6{u{x, t) — u)) 



-u{{uxx\u,h)P{u,h,t)}uu- (19) 



where {-{u, h) denotes the average conditional on a given -u, h. Therefore using the eq.(|T9D 
it follows that the P{h,u,t) satisfies the following, 



- Put= -iPj,^ - ^{u'P)f,,, - a{uPy^ - /c(0)P^^, + KMPuuu - ^{uPhj,+ 



+u{{u:rx\u,h)P(u,h,t)}uu- (20) 

This equation is exact for a given u and as it is clear the fingerprints of diffusion term is 
giving rise again to an unclosed equation for P{h,u,t). Accessing to the functional form of 
the conditional averaging {u^xlu, h) is one of the major difficulties in the formulation. From 
eq.(|l8D we see that P(u, h) = P{—u, h) which results to, 

{u\h) = 0, (21) 

for any given u. In fact for a given h the average of height gradient, u, is consequently 
zero. The identity proposed by eq.(l21|) is not restricted to any limiting asymptotic and 
is true in all regimes of the dynamical evolution of the surface. Also eq. (|20|) allows us to 
determine a dynamical equation for P{h — h). Doing so, we multiply the eq.(^) to u an 
integrate over u from — oo to +oo, from which we get, 

d^P{h,t) = ^{{-{{u') - {u'\h)))P{h,t)} 
oh 2. 

+^,{im - '^{n'\h))PCh,t)h (22) 

where h = h — h and the relation 7 = f (u^) is used. This is a Fokker-Planck (FP) equation, 
describing the time evolution of P{h,t). The drift coefficient in the FP equation is, 

= -|((n^) - {urh)), (23) 

and the diffusion coefficient reads, 

D^^^ =k{0)-u{u''\h). (24) 
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It is evident that aiming to obtain P{h — h), one should know the conditional average {u^\h). 
The equation has the following stationary solution: 

PstatXh) = exp{ £ dh'D('\h')/D('\h')} (25) 

Where N is the normalisation coefficient. Therefore to derive the moments of hight 
difference h — h i.e. {{h — /i)"), we need the conditional averaging (it^|/i). The simplified 
picture given by this equation is indicating that all the knowledge one would know in order 
to obtain the behaviour of PDF is buried in the functional form of one conditional average , 
i.e. {u^lh). Although simple but it is clear that the conditional average (u^lh) would have a 
non-trivial dependence on u and L in the limit of z/ — > 0. Instead of following this strategy 
however in the next section we follow another direct way of extracting the moments of height 
difference {h — h) in the strong coupling limit , i.e. i/ — > 0. 

III. THE JOINT CORRELATIONS OF HEIGHT DIFFERENCE AND HEIGHT 
GRADIENT BEFORE SHARP VALLEY FORMATION 

When a is finite, the very existence of the non-linear term in the KPZ equation leads to 
the development of the sharp valley singularities in a finite time and in the strong coupling 
limit (i/ ^0). In one dimension system is already in the strong coupling regime so starting 
from any finite value of i/ in large time system develops sharp valley singularities, Fig(2). 
Therefore one would distinguish between different time regimes before and after the sharp 
valley formation. Starting from a flat initial condition, i.e. h{x, 0) — 0, u{x, 0) = 0, which its 
evolution is given by inviscid KPZ equation, we know that after a finite time the derivative 
of function h{x,t) becomes singular. After this time scale the diffusion term is important, 
but we can neglect that before appearance of the singularities. So the equation governing 
the the evolution of generating function, Z{ii, A, t) before the creation of the sharp valleys 
is given by, 

Zt = i-f{t)XZ - iX^Z^^ - \^k{Q)Z + —Z^ - iaix{^)^ + fi^hMZ, (26) 

^ jJi jJi 
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where we have assumed the statistical homogeneity [Z^ = 0). Now we need the 7 which is 
given by 7 = ht. So invoking to the eq.(|l^) we get, 



l{t) = h = -{u'). (27) 

To evaluate (u^) we set A = in the eg. (p6|) , and find, 

Zt = ii^h,{0)Z, (28) 

which considering Z{fi, 0) = 1 as initial condition, its solution is, 

Z(/i,t) =exp(/i2A;,,(0)t). (29) 

On the other hand by definition we have (m^) = — ( ^ f^a'^^ )- So before creation of the 
singularities the second moment of height gradient behave as: 

{u^) = -2k,Mt, (30) 

so consequently, 

7(t) = -aA;,,(0)t. (31) 

Inserting eq.(PT|) in the equation (^) results to, 

d X d 

— A, t) = -iaX—Z{^, A, t) + ia-—Z{^, A, t) 

+ (/i'A;,.,.(0) - iah,{0)t\ - \^k{0))Z{fi, A, t). (32) 

We solve the eq.(|3^) with the initial condition Z(yU, A,0)=1, from which by expanding the 
generating function in powers of A and fi we can obtain the moments {{h — h)^) and 
{{h — JiYu^). We change the variable n to y = yU^, so the eq.(^) converts to the following 
equation, 

d d'^ d 

— A, t) = -2ia\y—Z{y, A, t) + ia\—Z{y, A, t) 
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+ {yhxiO) - laKMtX - X^k{0))Z{y, A, t). (33) 

Introducing the Fourier transform of Z{^, A, t) with respect to y as Q{q, A, t), it is simple to 
get the following evolution equation satisfied by the Fourier transform, 

—Q{q, A, t) = 2aXq^—Q{q, A, t) + 5aXqQ{q, A, t) - th^iO)—Q{q, A, t) + 
ot oq oq 

-iaKMtXQ{q, A, t) - X^k{0)Q{q, A, t), (34) 
with the initial condition 

Q{q,X,0) = ^ [e'y^dy = 5iq). (35) 
The eq. (0) is a first order partial differential equation which can be solved by the method 



of characteristics. The general solution of eq. ( pi] ) is written as 



( 1 2tfc,,(0)aA + y'-2«A;,,(0)aAtanh-i(gy^^^, 
^'2 KA^)aX 



J lOaAs + i tanh-i (s^g^) ^-2z/c,,(0)aA 

xexp|-l/2y^ V 2aA.2 - 1^.(0) 

+2mA;,,(0)tA - 2A2A;(0) - ^tanh-i(gyg^)y^-2zA;,,(0)aA | 

where 5^ is an arbitrary function of its arguments. Imposing the initial condition, given 
in eq . (p5|) , and introducing a; as. 



1 tanh-i(g^/-f^))7-2zA;,.,.(0)aA 
^ = o ^ — TT^r^ ' (37) 



2 fc,.(0)aA 

we obtain, 



g{X,uj) = 6{-lV2\r-^^tanh{J2ih,{0)aXto))exp[l/2 f 
2 V aX ^ Jo 



lOaXs + i tanh-^ k/^^l./-2iA;^^(0)aA - 2iujhJ0)aX - 2X'^k(Q) 

2aXs^-tk^M ^ ^' ^ ^ 

from which Q{q, X, t) is obtained as, 
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g(.,A,t)^.(A,t + .)expi-l/2/^ ^^J^^ 

2aW -iKM ^^^^ 
Inverse Fourier transforming of the solution in eq. (|39|) is straightforward, so after switch- 
ing to variable /i we get the following solution for ^(yU, A,t): 

A, t) = {l- tanh2(y'2iA;^^(0)aAt)) exp { - | ln(l - tanh^(y^2ifc^^(0)aAt)) 



5 I 

+ -tanh-^(tanh2(i/2zA;^^(0)aAt)) - A^A;(0)t 



16 l + tanh(y^2?A;,.,(0)«At) 2 V aA 'V ^ ^ J 

Since we are interested in moments {{h — h)^), so setting /i = in (|40|) and expanding 
the generating function in power of A we can obtain them all. For example expanding up to 
O(A^) it is easy to see that the first sixth order of moments are behaving as following, 

{{h - hf) = -h{h,{0)W - 6k{0)), (41) 
24 

{{h - hf) = --k^^iofaH', (42) 

{{h - hf) = -^kUO)Wt' - h'k^MWkiO) + It'kiOf, (43) 
1U5 D 2 

- 2288 4 

{{h - hf) = -{^k^Ofah'' + -kU^fah'km, (44) 

{[h - hf) = -^^^t\85783k^^{0faH^ + 299970k^^{0faH^k{0) 

+623700k^^{0faH^k{0f - 1247400/c(0)=^). (45) 

The important content of the exact form derived above is that through them the time 
scale of sharp valley formation can be found. Actually there is no guarantee that the follow- 
ing generating function can be derived from a physical probability density function. So one 
should first check the realisability condition, i.e. P{h — h,t) > 0. In fact the above moment 
relations indicate that different even order moments become negative in some distinct char- 
acteristic time scales. Closer looking in the even moment relations reveals that the higher 
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the moments are, the smaller their characteristic time scales become. So asymptotically the 
rate of decrease tends to for very large even moments, where t* = ( ^^^^''(o) )^^^ ( figure 
1). Therefore we conclude that after this time the far tails of the probability distribution 
function start to become negative, which is reminiscent of sharp valley creation. It means 
that after the characteristic time scale t* one should also consider the contribution of the 
relaxation term in the limit of vanishing diffusion in order to find a realisable probability 
density function of height field. In other words disregarding the diffusion term in the PDF 
equation is valid only up to the time scales in which the singularities are developed. Taking 
into account that a > 0, the odd order moments are positive in time scales before formation 
of sharp valleys. It means that the probability density P{h — h,t) in this time regime is 
negatively skewed. In figures [3-5] we have demonstrated the role of a on the time scale of 
creation of singularities. Substituting the fcxx(O) and k{0) in the expression of t* it gives us 
t* = (nY^^DQ^^a^'^^^a^^^ . Hence the smaller the a is the shorter the time scale of shock 
creation would be(see figures 3,4 and 5). 

IV. THE EQUATION OF JOINT PDF OF HEIGHT DIFFERENCE AND 
HEIGHT GRADIENT IN THE STATIONARY STATE 

Assuming a stationary state, we are interested in investigating the stationary solutions 
of eq. (pO|) in the limit z/ — > 0. Of course in the stationary state the sharp valleys are 
fully developed and one should also take care of the diffusion term in the PDF equation. 
The complicated term involved with the singularities , can be cured by using the method 
proposed in [73]. Let us define, 

G(u, h, t)= lim uluxxlu, h)P(u, h, t) 

= lirn ^{uxxix, t)5{h — h{x, t))6{u — u{x, t))), (46) 

,where in the last step in eq.(^6D we have used the definition of joint PDF P{u, h,t). As- 
suming spatial ergodicity, the average of the dissipative term can be expressed as, 
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T-'{uxx\u, h)P{u, h, t) = v{uxx{x, t)5{u — u{x, t))5{h — h{x, t)) 

= u lim — — / dxUxx{x,t)S{u — u{x,t))S{h — h{x,t)). (47) 

L~*oo L N J~L/2 

It is well-known that the w-field which satisfies the Burgers equation, gives rise to discon- 
tinuous or shock solutions in the limit z/ — > 0. Consequently for finite a the shock solutions 
are manifested in height field as a set of sharp valleys at the positions where the shocks 
are located, where they are continuously connected by some hill configurations Fig(5). It is 
noted that Uxx is zero at the positions where no sharp valley exists. Therefore in the limit 
^ only small intervals around the sharp valleys will contribute to the integral in the 
eq . (^71) ■ Within these intervals, boundary layer analysis can be used for obtaining accurate 
approximation of t), h{x, t). Generally boundary layer analysis deals with the problems 
in which perturbations are operative over very narrow regions across which the dependent 
variables undergo very rapid changes. These narrow regions (sharp valley layers) frequently 
adjoin the boundaries of the domain of interest, owing the fact that a small parameter {y 
in the present problem) multiplies the highest derivative. A powerful method for treating 
boundary layer problems is the method of matched asymptotic expansions. The basic idea 
underlying this method is that an approximate solution to a given problem is sought not as 
a single expansion in terms of a single scale, but as two or more separate expansions in terms 
of two or more scales each of which is valid in part of the domain. The scales are chosen, 
so that the expansion as a whole, covers the whole domain of interest and the domains 
of validity of neighbouring expansions overlap. In order to handle the rapid variations in 
the sharp valley layers, we define a suitable magnified or stretched scale and expand the 
functions in terms of it in the sharp valley regions. For this purpose, we split u and h into 
a sum of inner solution near the sharp valleys and an outer solution away from the sharp 
valleys, and use systematic matched asymptotics to construct uniform approximation of u 
and h. For the outer solution, we look for an approximation in the form of a series in i/, 

u= u°"* = uo + z/ui + 0(i^2). 
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^ /io + i,/^^ + O , (48) 
where uq and ho satisfy the Burgers and KPZ equation without dissipation term, 

UQt + au^UQ^ = -fx, 

hot - ^{dM' = /• (49) 

In order to deal with the inner solution around the sharp valley , let y = y{t) be the 
position of a sharp valley, define the stretched variable z — and let, 

u'^{x,t)=v{^^ + S,t), 

h"\x,t) = h{^^^ + 5,t). (50) 

The parameter 5 is a perturbation of the sharp valley position and v and h satisfy the 
following equations, 

pvt -a{u- i''q)v^ + avv^ = v^^ - f^{z, t), 

1^% - uuh^ + riu'^h^ - ^{hzf = uh^^ + i''^f{z, t). (51) 

^ ^ 2"" ' ^ ^ u+,U- are the height gradients in right hand and left 

hand sides of the sharp valley in the position y, ( see figure-2). We look for a solution in the 
form, 

V ^ Vq ^- VVi ^- 0{v'^), 

h^ho + uhi + 0{u'^). (52) 
To leading order we get for vq and hg, 

hoz = 0, 
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a{vo - m)woz = vo^^, (53) 

where we have assumed that variance of f{z,t) is a smooth function so that we can neglect 
its variation in the sharp valley region (/^ = 0). In other words we suppose that a » 
0{u){i.e.o- » the typical layer width). One can easily integrate the eq.(p3D and find that, 

ho = const, 

s asz 
vo = u - -tanh(^), 

in which s = s{t) = u^ — u^ is the shock strength. The boundary condition for this equation 
arises from the matching condition, 

lim = hmuT' = u±-. (54) 

Basically since h'^(P{z) = C — v p vl^{z')dz\ where C is the integration constant. So the 0(1) 
solutions of fp" give rise to 0(z/) solutions in /i™ field and only the integration constant is 
the 0(1) part of the solution of /iq". In fact the constant in nothing but the height value 
at the sharp valley position. Of course due to the height continuity at sharp valley position 
there is no boundary layer for KPZ equation meaning that the rapid changing term in the 
sharp valley layer occurs in hr^xx while the highest derivative in KPZ equation involves only 
hxx- The above analysis show that, to 0{v), eq.(|3) can be estimated as, 

\im v{{uxx\u, h)P{u,h,t)), 



hmz/ hm y^E / dzu^Au - u^'^ix^tMh - h'^'iy.^t)), 

J ■' 

hm / rfz<5(«-n-(z,t))5(/.-/i-(i/„t)), 

L->oo L iV ~^ J-oo 



liin Tm 22 / dzvoZS{u - vo)S{h - h{y,, t)), (55) 

3 
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where flj is a layer located at yj with width >> 0(z/). Using the eq.(p3D and, 

dzvQzz = dvQ-^ = advQ{vQ — u), (56) 

the 2;— integral can be evaluated exactly leading to the following result, 

i'{uxx\u,h)P{u,h,t) = a du dsQ{u,s,h,x,t) / dvo{vo — u)6{u — vq). (57) 

J J —00 •^"+f 

g{u, s, h,x,t) is defined such that g{u, s, h,t)dudsdhdx gives the average number of vallies 
in [x, X + dx) with u{y, t) G [u,u + du), s{y, t) G [s, s + ds) and h{y, t) G [h,h + dh), where 
?/ G [x, X + dx) is the sharp valley location. Eq. (p7|) indicates that the relaxation term in 
the strong coupling limit can be written in terms of some quantities which are defined in 
singularities (vallies). Indeed we characterise a sharp valley with four quantities, its location 
yj, its gradients at yj^^ (i.e. m+), yj^_ (i.e. u_) and its height from the h i.e. hj. Instead of 
u+ and «_ we have used the quantities u = ^++^- and s = s(t) = — u^. Later we will 
determine the time evolution equations which govern over these four quantities. 
Proceeding further we note that g{u, s, h, x, t) can be defined as, 

g{u,s,h,x,t) = {J25{u-u{yj,t)6{s - s{yj,t))6{h - h{yj,t))5{x - yj)). (58) 
j 

Due to statistical homogeneity the sharp valley's characteristics are independent of their 
location, so 

g{u,s,h,x,t) = pS{u,s,h,t)), (59) 

in which p = p{t) is the number density of shocks and S{u,s,h,t) is the PDF of 
{u{yo,t), s{yo,t), h{yQ,t)) conditional on ?/o being a shock location. Hence 

\im.u{uxx\u.,h)P{u,h,t) = —ap / ds I du{u — u)S{u, s, h,t). (60) 

l^-^O J-00 "'M+f 

Therefore the relaxation (dissipative) contribution in eq.(^) is written as, 

G{u, h,t) = —{ap / ds du{u — u)S{u, s, h,t))uu- (61) 
J— 00 •''^+3 
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So eg. (120|) is rewritten in the following form, 

- P^t = -7P^„ - - C,{UP\ - mPll^ + KMPuuu + G{u, h, t). (62) 

It is interesting that the the G-term comes from the relaxation term in the KPZ equation, 
but its explicit expression in terms of sharp valley's characteristics is proportional to a which 
is the coefficient of the nonlinear term in the KPZ equation. This indicates that without 
the nonlinear term in the KPZ equation there is no finite contribution for the diffusion term 
in the PDF equation when u ^ 0. Although this equation is exact for finite a however 
we can not solve it since the last term is not expressed in terms of P{u, h,t). Despite the 
existence of unclosed G-term still we can derive interesting information the moments using 
the above equation, we will study comprehensively the moments of height difference and 
height gradient, i.e. {{h — h)^{dxh)^) in the next section. 

But before it's worth of remarking that integration over h gives an equation for proba- 
bility density function (PDF) of u recovering the results in [73], 

Rt = -kxx{0)Ruu + {pa / / du{u - u)S{u, s,t)}„ (63) 

where R{u,t) = JP{u,h,t)dh and S{u,s,t) = J S{u, s, h,t)dh is the PDF of 
{u{yo,t), s{yo,t)), conditional on the property that yo is the shock (cusp) position. Be- 
cause of the statistical homogeneity, yo is a dummy variable. We finish the section with 
determining the relation between the density of vallies and the noise characteristics k{0) 
and kxx{0), that is, 

k..{0) = ^{s') (64) 

indicating that the forcing variance kxx{0) is related to the products of density of vallies p, 
(s^) and a. This relation has been found in [73] and its details are given in the appendix. 
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V. THE MOMENTS OF HEIGHT FLUCTUATION IN THE STATIONARY 

STATE 

Our goal is investigating the scaling behaviour of moments of height difference and 
height gradient in the stationary state. After sharp valley formation the lateral correlations 
produced by nonlinear term will grow with time. Dynamic scaling exponent z characterises 
the self-similar growth of the lateral growth. However in the stationary state the height field 
width saturates in the sense that lateral correlations are in average grown up to the system 
size. As it was explained, after the saturation the width scales as wo{L,t ^ L^) ~ L^. 
Having in our disposal the exact result x = ^ in one dimension [89], it would be natural to 
define P{h', u, t) as PDF of h', u, and t, where h' = and wq = L^^'^. Obviously P{h', u, t) 
is related to P{h,u,t) as P{h',u,t) = woP{h,u,t). From eq. (|62D it follows that P{h',u,t) 
in the stationary state satisfies the following equation, 

-iL~'^'Ph'u + -^L-^'^u^P),,^ - aL~'l\uP)y 

-k{0)L-^Ph'h'u + KMPuuu + G{u, h', t) = 0. (65) 

From eq.(]65|) it follows that the moments of {h'^u"^) satisfies the following equation in the 
stationary state, 

gn,m + KMrn{m - l)(m - 2){h"'u"'-^) + mn{n - l)k{0)L-\h"'-\"''-^) 

- mn-fL-^/''{h"'-^u'"'^) - n{m - 2)-L~^'\h"'-^u'^+^) = 0, (66) 

where 

gn,m = j dh' j duhVG{u, h', t). (67) 

Using the eq.(0), we can determine the explicit expression of gn,m in terms of the charac- 
teristics of vallies. Thus we find, 
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gn,m = —ctp / dh' / duh'^u"^{ ds I du{u — u)S{u, s,h' ,t))uu- (68) 
After integrating by parts it converts to, 

gn,m= —ctp / dh' / duh'"'m{m — l)u'^^'^{ ds / du{u — u)S{u, s,h' ,t)), 



= -apm{m - 1) T ds f dh'h"" f duvJ^-^{ T ' du{u - u)S{u, s, h', t)). (69) 

It can be integrated over u which leads to the following expression for gn^m (see appendix 
B), 

9n,m= ^{{hT{2u + sr-'[{m - l)s - 2u)]) 



+ {K{2u - s)'"-^[(m - l)s + 2m)])), (70) 

where h'^ = ^"J^ and is the height of given sharp valley . This means that the 
relaxation term in the strong coupling limit can be written in terms of only characteristics 
of the vallies i.e. u, s and 

At statistical steady state {t — > oo, (. . .)t = 0), and assuming the scale independence of 
gn,m^s, one can derive from eg. (|66[) , to leading order in the limit of L — * oo the following. 



^0 ^ ^ A;..(0)m(m-l)(m-2)' ^' ^ 

for m > 3. Putting n = then for instance results in the height gradient moments behaving 
as [73], 

ap 



2'"+3fc^^(0)(m + l)(m + 2)(m + 3) 



{{{2u - s)'"+2[(m + 2)s + 2u)]) - {{2u + sr+^[2u - (m + 2)s])}. (72) 
and for m = 3 we find, 
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where using the eq.(70) we have gn,3 = ^f{h''^s^). Assuming that gnm's are scale indepen- 



dent, at least in leading order, then one may argue that eq.(|75D builds up a relation between 



the n-th moments of height difference in term of the second moment Wo in a non-intermittent 
way. That is the n-th order moment is scaled linearly with order n. Rationalising the as- 
sumption of gn,m^ scale independence, one would look at the statistics of sharp valley en- 
vironment and the different processes involved in the sharp valley creation and annihilation 
which contribute dynamically. We will postpone the argumentation to next section and in- 
stead for the rest of this section we find out the consequences of scale independency of gn,mS. 
In the end of this section we will verify the requirements of positivity and normalisability of 
the joint P{h, u, t) in the stationary state. 

Eq.(|66|) also suggests that the amplitudes of height difference and height gradient mo- 
ments depend strongly on the singular structures in the theory, encoded in the functions 
gn,m (i.e. eq.(^). 

At this stage we find the finite size effects on the moments of Sn^m = {h'"'u"^). Defining 
e = as a perturbative parameter we find the structure functions Sn^m = {h''^u"^) 

perturbatively in terms of perturbative parameter e = -^172 as , 

Sn,m = {h "m™) = S^^Jj^ + €S^^^ + e^S'^^^ + ■ ■ (74) 



Invoking to eq.(|7^ and assuming the scale independency of gnm we get, 



K.mm + l)(m + 2)(m + 3)^'^''"+^' ^^^^ 

c.(i) ^ \ r 7n(m + 3) 

'"^ KMim + l)(m + 2)(m + 3) ^ fc,.,.(0)(m + 3)(m + 4)(m + 5)y"-i'"^+5 
an{m + 1) , 
^2fc,,(0)(m + 5)(m + 6)(m + 7)^"-i'-+7i' 



and. 



k^x[^)[m + l)[m + 2)[m + 



%{m + l)5?2i,^+4 - k{0)n{n - l)(m + 3)5i°22,m+2], (76) 
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and etc. For example the moments {{h — /i)") behave as, 

{(h - hr) = { \ n„ 3 



1 7n 2Q;n 



1 fc(0)n(n — 1) 7^n(n — 1) 



lla7n(n - 1) 40a^n(n - 1) n ,nrr-3/2n 



(77) 



Noting to the fact that the gn-1,5 and g-n-ij are not zero, therefore we conclude that the next 
to leading order correction for structure functions is 0(1/L^/^). Also eq. ( [77D shows that 
the amplitude of the correction terms to moments {{h — h)"') are related to the statistics 
of quantities which are defined on the singularities i.e. gn,m- Also it shows all of the 
moments {{h — h)^) (for even and odd n) exist and consequently the PDF of the h — h is 
not symmetric. However using the properties of the Burgers equation, it can be shown that 
only even moments of u are nonzero and all the odd moments vanish, hence the PDF of u 
is symmetric. 

Eq.([72l) enables us to determine the rate of growth of surface at the stationary state, i.e. ht. 
Using the KPZ equation one it is trivial to see that. 



lim ^{t) = ht = -{u') + up{s), (78) 

t—fOO Z 

where we have used the fact that (h^^) = —{ux) = — p(s) [73]. In the limit z/ — » the second 
term vanishes and, 

\im^{t) = ht = ^{u'). (79) 
Back to the eq. (|72[) , the ht is written in terms of properties of singularities as. 



„2 , 
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So in the stationary state, moments (m^s^) and (s^) would determine the growth rate. In 
other words for a given time in the steady state if one average the moments (u'^s^) and 
(s^), which is defined only on the vallies, can predict the rate of growth of surface. This 
provides a simple way to determine the ht in the stationary state. Now we prove that 
P{h',u,t — s> oo,L — s> oo) = P{h',u) is a positive and normalisable PDF. Regarding the 
positivity of the PDF we note that eg. ([BT|) indicates that P{h',u) satisfies the following 
equation in the limit of L — »• oo, 

rO nu-^ _ _ _ ^ 

kxxifi)Puuu = {ap j ds du{u-u)S{u,s,h,t))uu- (81) 

J-oo -^f+f 

Invoking to the method introduced in section 3, one may obtain, 

P{h',u) = -— — -- / ds du[— - (u-u) ]Soo{u,s,h'), (82) 

where Soo{u, s, h') is the PDF of vaUies with u, s and h' . Therefore positivity Soo{u, s, h') > 
implies that P{h',u) > 0. To check the normalisability of the P{h',u), consider the eq.([7TD 
with n = 0, m = 3, leading to, 

J dh'duP{h',u,t) = ^^j^^^Q^ 9o,3, (83) 

while explicit form of (70,3 can be found from the eq. ([70|) . The result gfQ.s = ^{^^) cooked 
up with the eq. ( P15|) gives, 

J dh'duP{h\ M, t ^ 00, L ^ 00) = 1 (84) 

, completing the proof of normalisability of the stationary state PDF P{h' , u). 

VI. STATISTICS FOR THE ENVIRONMENTS OF THE SINGULARITIES 

In this section we derive the PDF of quantities which characterise the sharp valleys. 
As is depicted in Fig. (2) and formerly is described, the evolution of the surface after the 
formation of singularities is determined by the dynamics of the sharp vallies and their sta- 
tistical properties. In a more quantitative sense one should attempt in characterising the 
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time evolution of h^^, u and s consequently. From there we show that to leading order of 



the expansion in terms of system size the gn,m^ do not depend on the scale L. Doing so we 
reach to such a level of describing the dynamics of the surface growth by which one may also 
trace the dynamics of the singularity enviornments. That makes a logical way to construct 
the pathway towards examining the statistical properties of singularity functions gn,m- 

The 

importance of such an analysis became clear in the last section where the determination 
of the finite size corrections to scaling in Sn^m — {h'^u"^) showed to be depended on the 
existence or lack of scale dependencies in the singularity functions gn,m- 
Let us turn to study the statistics for the environment of the singularities in KPZ equation. 
Defining ^{x, t) — —hxx{x, t) and let W{hy, u, s, x, t) be the PDF of 



conditional on i/j being a singularity position. In this section we will find the master 
equation governing the evolution of W{hv,u,s,^+,^-,x,t) in the limit of — > 0. Starting 
from the dynamical equation, 



K{x,yj,t) = -{h{yj + 2) + ^(%- - 2^^' 
u{x, yj, t) = ^iu{yj + |) + u{yj - |)), 



s{x, yj, t) = u{yj + 2 ) - ""(^i " 3 ^' 




(85) 



ut{z + x±) + auux±{z + x±) 



fx±{z + X±), 



(86) 



Ct{z + x±) + au{z + x±)Cx±{z + x±) + a'^C^iz + x±) 



fx±x± (^ H" X±), 



(87) 



we define, 



9{X+, A_, 77+, 77_, x+, X-, z, t) 



exp{—iX+h{z + x+) — iX-h{z + X-) — i^+u{z + x+) 



ili-u{z + X-) — ir]+^{z + a;+) — ir]-^{z + a;_)). 



(88) 



and. 
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then, 



3 



(89) 



pH^(/i+, X, t) = 

y dX+dX-dfi+dji-drj+dr]^ ^-ix+h+-iX-h--iii+u+-ii,-u--iv+e+-in-i- _ ^gg) 
We now derive equations for (0) and W. Using the equations (109-111) we obtain, 
(9),= -^A+((f< + /+)E^5(^-%■)) 

j 
j 

3 
3 

+ - V3¥ + E(^(^ - v3)Kt - tk)e) 

3 k 

-(ZiSiz-yi)5{t-U)e). (91) 
I 

S^{z) — ■^S{z) , the (yjk,ifc)'s are the points of singularity creations and the (yi^tiYs are 
the points of singularity annihilation due to collisions. Assuming homogeneity and using 
the following identity [73] , 

eS\z - yj) = {eS{z - yj)), - OxJiz - %■) - OxJiz - yj), (92) 

it follows that, 

(6),= -iX^{(\l + U)^e5{z-y,)) 

^ 3 

-^A_((|«^_ + /_)E^^(^-%)) 

3 
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- iiJ,+ {{-au+U:c+ + fx+)^OKz - Vj)) 

3 

-ilJ,-{{-au-U^_ + fx-)^Od{z - yj)) 

3 

- iTl+{{-ail - au+i^^ + f^^^^) e5{z - yj)) 

j 

3 

+Si - E2, (93) 

where Ei and E2 account respectively for singularity creation and collision events. These 
are given by, 

Si(A+, A_,/x+,/x_,r;+,r;_,a;+,a;_,^,t) = (^e5{z - yk)6{t - tk)), (94) 

k 

E2(A+, A_, 77+, ?7_,x+,x_,z,i) = (^e5{z-yi)5{t-ti)). (95) 

Invoking to the Novikov's theorem we have, 

(/±^) = (-iA±A;(0) - iX^k{x± - x^) - in^:kx{x± - x^) 

+ iV±kxx{O) + iVTix±-x^)){0), (96) 
{fx±0) = {-i>^±kx{x± - x^) - iii±kxx{0) + iiJ.^k^{x± - x^) 

+ iriTkxxx{x± - x^)){9), (97) 
{fx±x±0) = {-i>^±kxx{x± - x^) 

— iXzfkxx{x± — Xzf) — ilJ'^kxxx{x± ~ + iVi^xxxxi^) 

+ i'n^kxxxx{x± - x^)){9). (98) 
To average the convective terms we use 

iafj,±{u±^±Q) + ia'r]±{u±^±Q) — —a{u±Qx±) + iaX±{u'^Q) 

= -ia{e)x±^± + ia{e)r,± - iQ;A±(e)^±^±. (99) 
Notifying that. 
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Ox± = {i\±u± - iiJ,±Cx± - ir}±Cx±)0, (100) 
(49) = -(6)^^^^, (101) 
(e|e) = -(6),^,^, (102) 

Finally for (0)^ we find, 

(e)t=-yA+(e),,,,-^A_(e),_,_ 
-(e).,,_ - (e)._,, + -mv, + (e)._) 

-(A^A;(0) + AH(0) + 2A+A_A;(a;+ - x_))(e) 
-{lJ''ihx{0) + fJ-thxiO) + 2ij,+ii^k^^{x+ - x_))(©) 
-{vlkxxxxiO) + rf_k:cxxx{^) + 2ri+ri-k{x+ - X-)){e) 
-2(A+A(_ - A_//+)A;^(x+ - x_)(e) + 2(A+r7+ + A_r7_)A;a.a.(0)(e) 
+2(A+r/_ + A_r;+)A;^^(a;+ - x_)(e) - 2{r] 

V+l-''-)^xxx{x+ X-){Q) 

{aJ2u{yj: + e,_)5{z - yj)) + El - E2. (103) 

j 

for the term involving u{yj, t) we note that, 

u±{yj,t)9x± = {u{yj + x±,t)9)x± - ^{yj + x±,t)9 = i6'a;±A± - i6'^±, (104) 
li±(l/i, ^)6'a:=F = {u{yj + Xi, t)6')a:^ = i9^^X±, (105) 

thus, 

«(IZ^(%'^)(^^+ + ^^-)<^(^-%)) 
j 

j rv 

= Y((e>).+A, + (©)._A, + (e).,._ 
+(e)._A_-(e)^+-(e)^_). (loe) 

Combining the above expressions, on the subset A+ = A_ = I, a;+ = —X- — f,/xi = 
Ai+ + A*-, = ^"^2^" , (©) satisfies, 
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-Y(MO) + Kx)m - + 2^l)KM - 2(^ - iJ^DkU^m 

xxxx xxxx 

(0) + 2rj+r)-k 

xxxx 

+X{V+ + V-){kxx{0) + k,,{x)){Q) + 2fi2{v+ + V-)kxxx{x){e) + 
/xi(77_ - V+)kxxx{x){Q) + S2 - El. 

The El, E2 are evaluated at A+ = A_ = |, a;+ = — a;_ = |, /^i = /i+ + /X2 = . 
Changing to the variables {hy, u, s, ^_), we obtain the following equation for W, 

{pW{hy, u, s, X, t))t = ^u^pWh, + ^s'^pWh, - aspW^ 

+ ^^+pW + ^^.pW + api^+'W)^^ + ap{^JW)^_ 
+k{0)pWh,h. + k^^{0)pWm + 2p{k^^{0) - k^^{x))Wss 
+pkxxxxiO){W^+^^ + W^_^_) + 2pk^^^^{x)W^_^_ - 2k^{x)pWsh, 
-2k^^{0)pWh,^+ - 2k^^{0)pWh,^_ - 2k^^{x)pWh,^^ - 2k^^{x)pWh,^_ 
-2pt^^^.^{x){Ws^_ + Ws^+)t^^^{x)pWu^_ +k^^^{x)pWu^^ 

+Ci - C2- (107) 

The Ciihv , u, s, ^_|_, ^-,x, t) is defined such that, 

(i{hy.,u, s, ^-,x, t)dhy.dsdud^+d^_dzdt, (108) 

gives the average number of singularity creation points in [z, z + dz) x [t, t + dt) with, 

hy{x,yi,ti) e [K,hy + dhy), 
u{x,yi,ti) e [u,u + du), 

s{x, Di, ti) e [s,s + ds), 

e(z/i + |,^i) e [^+:^+ + d^+), 
e(i/i-|,ti) e [i.,i. + di_), 
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conditional on G {[z,z + dz) x [t,t + dt)) being a point of singularity creation 

(because of the statistical homogeneity, z is a. dummy variable). (2{hvj,u,s,C,+,C,-.,x,t) is 
defined such that, 

C^ihvjjU, s, C,-,x, t)dhy.dsdudS,+d^^dzdt, (109) 

gives the average number of singularity collision points in [z, z + dz) x [t, t + dt) with, 

hv{x, 1/2, h) e [K, + dh^), 
u{x,y2,t2) G [u,u + du), 
s{x,y2,t2) e [s,s + ds), 

e(Z/2 + |,t2) G [e+,^+ + rfe+), 
^(l/2-^,t2)G[^,^+rfe-), 

conditional on {1/2, t2) G ([-2,-2 + (i-z) x [t, t + dt)) being a point of singularity collision. 
Now we rescale as h'^ = jtj2-, so the equation ( |1 07] ) changes to, 

+ap(e+^iy')€+ + «P(e-'W^')?- + ^-'MO)pW-/:^,i + kU^)pWU 
+2p(A;,,(0) - A;,..(x))iy;, + pk^^U^){W'^^^^ + 
+2pA;,,,,(x)iyg'_^_ - 2L-^'^k^{x)pW',^, - 2^,(0) L''/^ pW^,^^ 
-2hML^'^'pWk^_ - 2kUx)L-"^pWi.^^ - 2kUx)L-^'^ pWi,^_ 
-2pk^^^{x){W'^^_ + W'^^^) - k^^^{x)pWL^_ + k^^^{x)pW'^^^ 

+C'i - C2 = 0. (110) 

In the limit of large L or L — 00 the leading terms are, 

-aspW^ + -^^pW + -^.pW 
+api^^'W%^ + api^JW')^_ + k^MpWLn 
+2p(A;,,(0) - k,,{x))W^, + pk,,,,{0)iW^^^^ + W^^J 
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hMpW^^_ + k,,,{x)pW;,^^ +([-('2 = 0- (111) 

To find the gn,m we multiply the above equation by K'^vJ^s^ and integrating over 
/i', s, ^+ and ^_ we have, 

+2p(p - l)p[kU^) - k,,{x)]{}{:irsP-'') 

+m(m - l)kUmK^n-^-'s^) + - Qi^^ = 0, (112) 



where, 



Q^nlp = / l^'^lTs^'^dh'dudsdi+di-. 



Using the following identity. 



we find. 



—h'^{x,yj,t)^-^, (113) 

^six,y„t) = ^{^+ + ^_), (114) 
d s 

— M(a;,y^-,t) = -, (115) 



So in the limit of L — > 00 we have 
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+2p{p - l)p[hM - hAx)]{h':u'^sP-') 

+m(m - l)KMp{h':u^-'sn + - Q<^1^ = 0. (117) 

Assuming a stationary solution for the dynamical equation ( |107| ) governed over W and 
rescaling hy as h[ = jtj2 in the resulting differential equation we reach to eq. (|110|). Of course 



the mentioned equation is dependent on scale L but being interested in the limit of L — oo 
results to eq. (|111|) which is free of the explicit scale dependent terms in the leading order 



. However we are faced with two very complicated terms namely (i and (2 which would 
be analysed . The origin of these terms are related to processes of sharp valley creation 
and annihilation. We argue that these processes are basically involving local interaction 
between nearby sharp valleys and effects of forcing, which its spatial correlation is assumed 
to be much less than system size, so they essentially would not carry any information about 
system size. In this sense eq. ( P-ll| ) encodes the fact that the probability distribution is a 
scale invariant function of its argument h'^ = -jjj2 in the leading order. The above property 



is deciphered in eq.( 117 ) too but this time it is translated in terms of the scale independence 

of 5'n,mS. 

Also the equation for W enables us to find the time evolution of the sharp valley character- 
istics. For example multiplying the equation ( |107| ) by and integrating over all variables 
we can derive the increasing rate of mean height of the singularities and noting that, 

d s 

— hy{x,yj,t) = --, (118) 
^s(x,i/„t) = ^(e+ + e-), (119) 

we get, 

j^{K){t) = -^{Au'-s'). (120) 
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VII. CONCLUSION 



We study the problem of non-equilibrium surface growth described by forced KPZ equa- 
tion in 1-1-1 dimensions. The forcing is a white in time Gauaaian noise but with a Gaussian 
correlation in space. Modelling a short range correlated noise we restrict our study to the 
case when the correlation length of the forcing is much smaller than the system size. In the 
non-stationary regime when the sharp valley structures are not yet developed we find an 
exact form for the generating function of the joint fluctuations of height and height gradient. 
We determine the time scale of the sharp valley formation and the exact functional form 
of the time dependence in the height difference moments at any given order. Investigating 
the stationary state we give a general expression of the mixed correlations of height and 
height-gradient at any order, in terms of the quantities which characterise the sharp valley 
singular structures. Through a careful analysis being done over the behaviour of the sharp 
valley environment, wc decipher the general finite size corrections to the scaling of an arbi- 
trary nth moment, i.e. {{h — h)"'), at any order. Recently Marinari etal. [23] have obtained 
the corrections to the leading order scaling in dimensions D — 2,3,4, in a high resolution 
simulation on the RSOS discrete model which is beleived to be in the universality class of 
the KPZ equation stirred with a white in time Gaussian noise and delta correlated in space. 
Hence they get, 

Wn(L) ~ A„L">^(1 + S„L-). (121) 

Irrespective of the dimension and moment order n, they observe the same sub-leading 
exponent uj always very close to unity (see also [91,92]). Through our calculations wc succeed 
to obtain the finite size corrections analytically. However we have to remark that, due to 
working with finite correlated forcing a firm comparison between our results and numerical 
simulations is not possible. More precisely, in the present paper the limiting of z/ — is 
taken into account only when a is finite. Still the forcing correlation length is much smaller 
than the system size and height correlation length. But the limiting of cr — > is a singular 
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limit in our calculations, and moreover, it is not a priori clear that the limits of — > 
and cr ^ commute at all. However due to the scale independence of gn,m's, the eq.([77D 
shows the general correction terms for nth order moment, all having the same sub- leading 



exponent uj = 1/2. The amplitudes An and Bn in eq.(121) are given explicitly in terms 
of the functions gn,m defined on the sharp valley singularities. The next step, left for the 
future, would be the calculation of gn,m^s in terms of few known parameters, i.e. the forcing 
and diffusion coefficients. 

Our analysis enables us to find the stochastic equations which are governed over the dy- 
namics of quantities characterising the cusp singularities too. This translates the stationary 
non-equilibrium dynamics of the surface in terms of the dynamics of singularities in the sta- 
tionary state. When the system crosses over the time t*, after which the first singularties are 
formed, it would be an important study to analyse the shape deformation of non-stationary 
height PDF P{h',t) in time. We believe that the analysis followed in this paper is quite 
suitable for the zero temperature limit in the problem of directed polymer in the random 
potential with short range correlations [88]. The same method applied to KPZ equation in 
higher dimensions would be definitely one of the consequent goals of the present work. The 
main message which might be encoded in the present work is the importance of the statis- 
tical properties of the geometrical singular structures for understanding the strong coupling 
regime of Kardar-Parisi-Zhang equation in higher dimensions. 
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APPENDIX A: AN ALTERNATIVE METHOD FOR DETERMINING OF THE 
MOMENTS OF HEIGHT FLUCTUATION BEFORE THE FORMATION OF THE 

SINGULARITIES 



In this appendix we give the details of calculations of the scaling behaviour of moments of 
height difference before the formation of singularities. We know that the generating function 
X,t) satisfies the following equation when (z/ — > 0), 

Z, = t^XZ - ^Z^,. - X^kiO)Z + —Z^ + fi^KMZ. (Al) 

Let us write Z{fi, X,t) as follow, 

Z{fi, X, t) := (1 + A{t)X^ + C{t)X^ + F{t)X^ + G{t)X^ + J{t)X^ + M{t)X^ 
+B{t)X^? + D(t)AV' + E{t)X^i^ + iJ(t)AV + K{t)X^i2'' + L{t)X*iJ^ + 
N{t)X^fi^ + P(t)AV + exp {{~X^k{0) + fi^KM)t). (A2) 

Now expanding Z{^, X, t) as a series of /x, A and substituting it in the equation eq. (PB]) , we 
equate the terms in different orders of /i, A ending with some coupled differential equations 
governed over the coefficients introduced in the definition of Z in eq. (PUD . So we have, 

^^A{t)=%aB{t), (A3) 

^S(t) = -2iak^M^t\ (A4) 

^C{t)=iaD{t), (A5) 

^D(t) = -Aiak,MtB{t) - 2taE{t), (A6) 

= 0, (A7) 

^F{t)=iaH{t), (A8) 

^^H{t) = -Aiak^^{0)tD{t) - 2iaK{t) - 2iak^^{0)H^ A{t) , (A9) 
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^K(t) = -2taKMH^B{t) - 8taKMtE{t) - 9taQ{t), (AlO) 

= 0, (All) 

By solving these differential equations with the 

initial conditions that A{t), B{t),C{t), D{t), E{t), F{t), H{t), K{t),Q{t) are zero at t = 0, 
we find, 



A{t) = -a^KMH\ 




(A12) 


B{t) = -'^taK.iOft^ 


(A13) 
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(A14) 


D{t) = Aa2fc.,(0)¥, 


(A15) 


E{t) = 0, 


(A16) 


F{t) = aXj'. 
V J 2520 "^"^ ' 

Hit) = ^taXj', 


(A17) 
(A18) 


K{t) = -la'kU^)H\ 


(A19) 


Qit) = 0. 


(A20) 



By replacing these expressions in eq.(|A^) we find Z(n,X,t) as a function of n,X,t ex- 
plicitly without any unknown terms or expressions. Now if we expand the original form of 
generating function Z{fi, X, t) as a series in /x, A we find, 

Z{i^,X,t) = {expiX{{h-h))+ii2{d^{h-h))) = 



720 ^ 120' ' ^ 48' ' ^ 36 

-^(h - /i) V/i'A^ - —(h - hfu^iX'' - ^{h - JifX"" - — zuV 
48^ ' ^ 120^ ' ^ 720^ ' 120 

- h)u^^/X - —i{h - /i) V/i^A' - —i{h - /i) V/i^A^ - —i(h - /i) Va^ 
2^ X2 12 24 

- /i)^A^ + -t-wV + - /i)wVA + 7(/i - h)\'fJ,^X' 
1/0 24 o 4 

_ h)\fix^ + ^{h - ^^x" + i^wV + - h)u'^?x 

b 24 b 2 
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+ -i{h - Kfu^X^ + -i{h - hf\^ - -n^ - {h - h)ufiX 
2 6 2 

~{h-hf\^-iun-i{h-h)\ + l. (A21) 



Equating the coefficients of eg. (|A2|) and eg. ( |A21| ) proportional to the same powers in fj, 

and A and replacing the expressions of A{t), B{t),C{t), D{t), E{t), F{t), H{t), K{t),Q{t) it 
results to the same expressions given before, i.e. 

{{h - h)') = -h{K^{0)W - em), (A22) 
24 

{{h-hf) = ~-KMW, (A23) 

{{h - ly) = -^k,,{0)^aH' - ^t'k,40)'a'k{0) + h'k{Of. (A24) 



APPENDIX B: THE PROOF OF RELATION BETWEEN p AND {S^) 



We consider the statistical steady state i.e. Rt = so that the eg. (|63D can be written as 
follows, 

1 ap f'^ /""^f 

Ruu= , (n\ ^^^'^^ ^ ir (n\ ^l I ^ du{u - u)S{u,s,t))u. (Bl) 

We integrate the eg. (|B1|) w.r.t u and find 



oo Ju 



ap /"O /"""f 

Ru = -, — 77-r / ds I du{u — u)S{u, s,t). (B2) 



00 Ju 



kxx{^) 

At the large time limit {t 00) we denote R and S as R^o and Soo{u, s). So 



ap 







Roo = -, — 77^/ du ds du{u-u)Sooiu,s). (B3) 

'^xx{'^) J -00 J -00 "'«+f 

To determine the Roo we define function K{u) as the following. 
Differentiation the above eguation with respect to u gives us. 



^K{u) = / / ^ du{u-u)Soo{u,s) 

au kxx[^) J -00 Ju+^ 
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ap 



Now we integrate the eg. (p5|) over u from — oo to u and find, 



/ du—K{u) = - — ^ / du I ds [ ^ du{u — u)Sooiu, s 

J-oo du krriO) J-oo J -oo Ju+S. 



du / cts — Ooo(,M ,s) ; — -— / du ds—boo[u-\ — ,S). (Bo) 



—oo J— CO 4 2 2/c^^(0) J —oo J —oo 4 2 

Then we will find, 

K{u) - K{-oo)= Roo{u) + — — — - / ds du—Soo{u,s) 

^"'XX^Uj J —oo J-oo 4 

'^'^ ^ ds f ^ du — Soo{u,s). (B7) 



2k^^{0) J-oo J-oo 4 

According to the definition of K{u) we see that K{—oo) ( the shock probability density 
function goes to zero in this limit) and therefore we find the following relation between K{u) 



and Rooiu), 



K{u) = R^{u) + ds r " dU^S^iu, s). (B8) 

2fcxa;(,U) J-oo 4 

Using the eqs.(p4|) and (|B8|) we find explicit relation between the i?oo and Soo{u, s) as follows. 



^-(") = -^71^ f T"' ^^[t - - ^)']Soo{u, s). (B9) 
Assuming 5*00(^^,3) > , it becomes evident that above integral would give a realisable 



portability density for height gradient, that is R^o > 0. For finite a the eg. (p9|) gives us 
the PDF of height gradient in the KPZ eguation in the strong coupling limit. The function 
Roo{u)) enables us to determine the relation between the vallies density p and k^^{0). We 
would integrate over u from R^o so we define another function Ki{u) such that, 

KAu) = r ds du\-u - ^^^^—^]S^{u, s), (BIO) 

^ ' 2/c^^(0)i-oo A+f U 3 J oov , y ) 
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where differentiation Ki{u) with respect to u gives, 



ap 



2/^x0; (0) 

Now integrating the above equation over u from — oo to +00 gives, 



/+00 
duRoo {u 
-00 



2A)a;x(0) 



/+00 rO ^2 ^ ^3 

/+00 rO g 

-00 d^J_j4ji^-2^ + ^]S^i^ 



2fca;a;(0) 

Using the fact that Ki{+oo) — Ki{—oo) — 0, we obtain. 



in which the sum of the terms in right hand side gives, 

f + CXD 



^ 12A;,,(0) 

Thus from the requirement that be normahsed to unity, we get. 
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APPENDIX C: DERIVATION OF THE FINITE CONTRIBUTION OF 
RELAXATION TERM IN THE STATIONARY STATE 



In this appendix we give the detail of calculations of gn,m in the equation eq. (|70D . For 
computing the gnm we introduce, 

rO POO pu+§ n,m-l m 

K(u) = -m{m-l)ap ds dh!h"^ du{ u }S{u,s,h\t). (CI) 

J-oo J~oo "'«-f m — 1 m 

By differentiating K{u) and integrating in the whole range of u we have, 

dK{u) 



dK[u) , 
/ — 1 du = gr, 



ds dh'h"" du{ {u- -)}S{u- -,s,h',t) 

-oo J-oo J-oo m m — 1 2 2 



m— 1 



/U fOO fOO y"^ II S S 

ds dh'h"" du{ {u + -)}S{u + -,s,h',t). (C2) 
-oo J —oo J —oo Tfl Tfl — 1 2 2 

Since K{+oo) = K{—oo) = 0, left hand side vanishes, so, 

pO poo ^ ^ POO [y^ _|_ _|_ s-jm-l ^ 

gnm = rn{m — l)ap j ds j dh'h'^ I du[ — — u\S{u,s,h\t) 

J —oo Tfl Tfl — 1 



oo J —oo 



pO poo poo (y — (y — s_\m-l 

-m{m-l)apj ds j dh'h"" du[- ^- ^ u]S{u,s,h\t) 

J —oo Tfl Tfl — 1 



oo J-oo 



n pO rOO POO 

/ ds du dh'h"'{{2u + s)""-^ 

J —oo J —oo J —oo 



{[{Tfl - l)s - 2u] + (2m - s)™-^[(m - l)s + 2u]}S{u, s, h', t), (C3) 
which finally leads to eq. (|70D . 

APPENDIX D: DYNAMICS OF QUANTITIES WHICH ARE DEFINED ON THE 

SHARP VALLEYS 

In this section we determine the equation of motion for u{yj,t), s{yj,t), h{yj,t) along 
the sharp valley which is located at position yj at time t. Using the KPZ equation and its 
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differentiation by x around the sharp valley at position y^ , one can find a set of equation for 
h+iVj, t) = lime-.o+ h{yj + €, t), h^{yj, t) = lim,_o+ h{yj-e, t), u+{yj, t) = lim,_o+ u{yj+€, t), 
and U-{yj, t) = lim£^o+ u{yj — e, t) as follows, 

Ktiyj,t) = ^uliyj,t) + fiyj,t), (di) 

h^,{yj,t) = V + (D2) 

, t) = -au+ivj, t)u+^{yj, t) - f^{yj, t), (D3) 
U-tiVj, t) = -OiU^iVj, t)u-xiyj^ - f^{yj, t). (D4) 

To determine the -^{u, s, h} we use the following identity [73,84], 

ju+{yj, t) = ^u+^{yj, t) + u^tiyj, t) 

= au{yj, t)u+^{yj, t) - au+^{yj, t)u+{yj, t) - f^{yj, t) 
ot 

= - 2 ^) - fxiyj^ t), (D5) 

where u = Similarly, 



d _, . a 
—UlVi, t) = — 

dt ' 2 



= n^yj^ t)u-Ayj, t) - f^iyj, t). (D6) 



These equations can be re-written as, 

t) = -JS{U+^ - M_ J - 

where will give the equations for u and s. Since = —h^ we write the above equations in 
term of curvature of the surface in the right and left sides of the sharp valley at position yj 
as, 



d _ a 



LX 

^(yj^'t) = -jSih+a^x - h-xx) - fx, 



d Ot 

j^s{yj:t) = -s{h+,, + h_,,). (D8) 
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For determining the time evolution of the h — h — h, we use the KPZ equation by which 
one can easily show that h+{yj,t) and h-{yj,t) satisfy, 

jh-{yj.t) - + (D9) 

By definition we have j^h{yj,t) — + j^h-), so using the equation for h+ and 

and 

|%-,i) = -^(4tZ2-.^) + /-7, (Dll) 
where h{yj, t) — h{yj, t) —h and ht — 7. 

Therefore in summary we have the following set of equations for given sharp valley in 
the KPZ problem in the limit v ^ 0, 

dyj 

—— = au, 
dt 

dt Oi 

j^h{yj,t)^--{Au'-s')+f-r (D12) 
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Figure Captions 



Figure. 1 Behaviour of the time scales in which the moments {{h — /i)^") become nega- 
tive in terms of n. The square points are calculated according to eg . (^0]) while the solid line 
is the fitting curve asymptotically tending to \{ a2k^\o) Y^^ ■ 

Figure. 2 In the upper graph the sharp valley solutions in KPZ equation are demon- 
strated while in lower one the corresponding shock structures in Burgers equation are 
sketched . The variables characterising the cusp, namely h^- , h^j^ and h are shown. 

Figure. 3 Different time snapshots of gradient configuration whithin system size, i.e. 
—dxh verses x. The time scale for shock creation is demonstrated for a L. The solid 
points are showing the jumps in the height gradient. J. Bee [93] . 

Figure. 4 Different time snapshots of gradient configuration whithin system size, i.e. 
—dxh verses x. The time scale for shock creation is demonstrated for a < L. The solid 
points are showing the jumps in the height gradient. J. Bee [93] . 

Figure. 5 Different time snapshots of gradient configuration whithin system size, i.e. 
—dxh verses x. The time scale for shock creation is demonstrated for a <^ L. The solid 
points are showing the jumps in the height gradient. J. Bee [93]. 
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